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Background:  Delay  Transition  Using  Non- 
Equilibrium  C02  to  Suppress  the  Second  Mode 


PROBLEM:  In  hypersonic  flight,  heating  loads  are  typically  a  dominant  design  factor 

Turbulent  heat  transfer  rates  can  be  about  an  order  of  magnitude  higher  than  laminar  rates  at 

hypersonic  Mach  numbers 

A  reduction  in  heating  loads  by  keeping  the  boundary  layer  laminar  longer  means  less  thermal 
protection  needed  and  hence  less  weight  to  carry,  or  conversely  more  payload  deliverable  for  a 

given  thrust. 

OBJECTIVE:  Delay  transition  from  laminar  to  turbulent  flow  in  the  boundary 
layer  of  a  slender  hypersonic  body  by  using  nonequilibrium  C02 


Transition  in  high  Mach  numbers 
occurs  through  the  Mack  mode  - 
amplification  of  acoustic  waves 
traveling  in  the  boundary  layer 


Molecular  vibration  and 
dissociation  damp 
acoustic  waves 


e? 


At  relevant  conditions,  C02 
absorbs  most  energy  at  the 
frequencies  most  strongly 
amplified  by  2nd  (Mack)  mode 


Inject  C02  to  delay  transition  in  air  flows  of  interest 
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Amplification  rate/cycle  or  absorption  rate/wavelength 
cOj/f  or  ar_.,,,AA 


Background 


*  For  C02  the  broad  sound  absorption  curve  peak  coincides  with  the  amplification  peaks 

*  This  coincidence  is  most  pronounced  at  enthalpies  of  ~10  MJ/kg 


6- 


6 


9 


O- — 0--0 


Bending 
modes  (OvO) 


Asymmetric 
modes  (00d) 


O  =  3380 
K 


O  =  960  K 


O  = 1993  K 


Typical  T5 
T*  *  1700  K  (C02) 


Aib  “  1  MS 

(Camac  1966) 

Plots  from:  Fujii,  K.,  Hornung,  H.G,  "Experimental 
Investigation  of  High-Enthalpy  Effects  on 
Attachment-Line  Boundary  Layer  Transition/'  AIAA 
Journal,  Vol.  41,  No.  7,  July  2003 
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Baseline  configuration 


Free-stream  parameters  correspond  to  Run  2540*  in  GALCIT  T5  shock  tunnel 


M  =  5.3 

DO 

r*  =  1323.77  K 

oc 

p*  =21993  Pa 

*  DO 


p =  0.05788  kg/m3 
U'  =  3866  m/s 
T'  =293  K.  T  IT  =  0.22 

w  7  tn  *  'DC 


u*  =  4.897  x  10~5  Pa  •  s 

r  DO 

L*  =1  m 


Gas  is  injected  with  the  total  mass  flow  rate  ranging  from  3  g/s  to  13.5  g/s 


*Parameters  are  determined  using  Mm,  Tx,  and  p^  reported  by  Wagnild,  R.M.  et  al.  (AIAA-20 10-1 244) 
and  perfect-gas  model  with  Pr=0.72  and  y=1.4 
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Numerical  approach  for  mean  flow 


In-house  Navier-Stokes  code  HSFlow* 

•  Perfect  gas  of  Pr=0.72,  y=1 .4 

•  Sutherland  viscosity-temperature  dependence 

•  Implicit  second-order  finite-volume  method 

•  Shock-capturing  scheme 

•  Third-order  WENO  for  advection  terms 


597x649  grid  with 

•  50%  clustering  in  the  boundary  layer 

•  Clustering  near  the  injector 


*Egorov,  I.V.  et  al.,  Theor.  Comput.  Fluid  Dyn.,  Vol.  20,  No.  1,  2006,  pp.  41-54. 
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Stability  analysis 


Local-parallel  stability  computations 

•  Third-order  Rungge-Kutta  scheme  for  integration  of  stability 
equation 

•  Gramm-Schmidt  orthogonalization  procedure 

•  Eigenvalues  are  calculated  using  a  shooting/Newton-Raphson 
procedure 

Disturbance-  q(t/)exp(ia£  +  i/3z  -  iut) 


For  temporal  problem x)  is  complex,  growth  rate=u;. 

For  spatial  problem  a(u>,  [3,  ^)is  complex,  growth  rate=cr  —  —a. 

X 

N-factors  N(x,  u;,  j3)  —  a  (u;,  /?,  x)dx 


*Egorov,  I.V.  et  al.,  Theor.  Comput.  Fluid  Dyn.,  Vol.  20,  No.  1,  2006,  pp.  41-54. 
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Boundary  layer  thickness 


baseline  configuration  with  injection  rate=1 3.5  g/s 
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Mean-flow  profiles 


U-velocity 

0.010  -TT - n - n — . — 1-| - rn - rn - m - rn - m - rn - 


0.008  - 


0.006  - 


0.004  - 


temperature 


Baseline  configuration  with  injection  rate=13.5  g/s 
•Thick  region  of  cold  dead  flow  near  wall 
•Slow  relaxation  downstream 
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We  are  dealing  with  acoustic  instability 


Mode  3=Mack  fifth  mod 
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AFFM^j  Temporal  instability  in  the  relaxation  region  for 
. baseline  configuration 

-  - (injection  rate  13.5  g/s) 


x=0.3 


•There  are  seven  unstable  modes! 

•Mode  2  (Mack  fourth  mode)  has  maximal  local  increment 
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We  are  dealing  with  acoustic  instability  (cont’d) 


Moae  b,  x=u.zyo,  a=on  i 


Dispersion  relation  from  WKB  analysis' 


kA-ZVA  ■ 


i 


\{ptU  —  Uj)2  M2  7T 

- - - -  —  1  ay  —  — b  7r  m 

T  4 


rn  =  0.  1. ... 


Acoustic  modes  are  formed  in  the  waveguide  between  the 
wall  (y=0)  and  the  relative  sonic  line  y=ya:  U(ya)=cr-a(ya) 

‘Guschin,  V.R.,  &  Fedorov,  A.V.,  Fluid  Dynamics,  Vol.  24,  No.1, 1989 
“Guschin,  V.R.,  &  Fedorov,  A.V.,  NASA-TT-20683,  April  1990 
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Injection  affects  N-factors  of  Mack 

second  mode 


AM<PTH 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


13 


N  Factor 


Perfect  gas  model  captures  basic 

trends 
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x,  m 


N-factors  including  real  gas  effects  (Wagnild  et  al.  AIAA-2010-1244) 
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How  to  improve  the  injector 
performance? 


*  Negative  slope  may  compensate  the  blowing  effect  on 
the  displacement  thickness 


*  Injection  of  zero  total  mass  addition  may  help  to 
reduce  the  relaxation  region 


v  ▼  ▼  Blowing 

Suction  0f  co2 

of  air 
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Injector  of  conical  shape 


/\MO>TH 


5-deg  half-angle  sharp  cone  with  the  injector  having 
the  slope  0=0C-A0 


A r  «  6 [tan  9  —  tan(0_  —  A0)]  «  bA9 
A  r  A  r 


A  b  = 


tan(0  A(9)  tail  9 
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Mean-flow  pressure  for  conical  injectors 

(injection  rate  13.5  g/s) 


004 


003 


003 


Injection 

region 


17 


AM<PTH 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


afb^1  Mean-flow  boundary  layer  thickness  for 

conical  injectors 

(injection  rate  13.5  g/s) 


Boundary-layer 
thickness  increases 
with  A6  in  the 
relaxation  region 
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shaping  A0=1°  without  injection 


Effect  of  shaping  is  local  and  small 
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Shaping  A0=1°  with  injection  of  13.5  g/s 


x,  m  x,  m 

Conical  injector  with  A0=1° 

•  Slightly  destabilizes  flow  in  the  near-field  relaxation  region 

•  Slightly  stabilizes  flow  in  the  far-field  relaxation  region 

/\MO>TH 


Run  2540 


13.5  g/s 
13.5  g/s,  1  deg 
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Shaping  A@=1°  and  injection  6.75  g/s 


x,  m  x,  m 


Shaping  of  A6=1°  produces  small  effect  on  stability  of 
flow  at  the  injection  rate  6.75  g/s 
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Shaping  of  A0=3°  without  injection 


Effect  of  shaping  is  local  and  small 
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Summary  plot  of  N-factors  for  conical 

injectors  (injection  rate  13.5  g/s) 


Shaping  with  A0=1°  and  3° 

•  Slightly  destabilizes  flow  in  the  near-field  relaxation  region 

•  Stabilizes  flow  in  the  far-field  relaxation  region 
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Injectors  of  cylindrical  shape 
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Mean  flow  for  cylindrical  injector 

(injection  rate  13.5  g/s) 


Injection 

region 


Axial  velocity  near  injector 
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Mean  flow  for  cylindrical  injector 

(various  injection  rates) 


Wall  pressure  Boundary  layer  thickness 
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N-factors  for  cylindrical  injector 


This  shaping 

•  Slightly  destabilizes  flow  in  the  near-field  relaxation  region 

•  Stabilizes  flow  in  the  mid-  and  far-field  relaxation  region 
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Summary  plot  of  N-factors  for  cylindrical  injector 

(various  injection  rates) 
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Jgy>  Estimates  of  the  injection  effect  on  the 
=  transition  onset 


m/m0 

Cylindrical  injector,  m0=13.5  g/s 
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Suction-blowing  of  zero  mass  addition 


yi 

J 

^rflTh' 

-//2xn  bio,,ina  1/2  x' 

—l/2<x  —  x()<l/2 
Z  —  (x  —  X0)27T  /  l ,  —  7T  <  Z  <  7T 


q(z)  = 


l  . 

qn  —  sm  z 
0  2tt 

l 


m 


+ 


=  9»9_/ sin  zdz  =  %l  /  * 


m  — 


o 

o 


—  f  sin 

0  9tt  J 


-90^  /  *- 
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Mean  flow  for  suction-blowing  system 

(m+=6.75  g/s) 


Injection 
region 
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Wall  pressure  distribution  for  suction-blowing 

system  (m+=6.75  g/s) 
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N-factors  for  suction-blowing  system 

(m+=6.75  g/s) 
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Conclusions 


•  Injection  induces  a  cold  dead- flow  layer  in  the  downstream 
relaxation  region 

•  The  near- wall  flow  behaves  as  a  wave  guide  which  can  support 
several  unstable  modes  of  acoustic  type 

—  The  most  unstable  is  the  Mack  second  mode 

—  The  phase  speeds  of  instability  are  close  to  those  of  slow  acoustic  waves 
in  the  free-stream 

—  Instability  frequencies  are  several  times  smaller  than  in  the  no  injection 
case 

—  This  may  lead  to  dramatic  increase  of  receptivity  to  free-stream  noise 
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Conclusions  (cont’d) 


•  The  eN  computations  for  baseline  configuration  showed 

—  Injection  leads  to  destabilization  of  the  near- field  region,  stabilization  of 
the  mid-field  relaxation  region,  and  destabilization  of  the  far-field 
relaxation  region 

—  The  level  of  these  effects  essentially  depend  on  the  injected  mass  flow 
rate 

•  The  injector  shaping  considered 

—  Does  not  stabilize  the  near-field  flow  at  sufficiently  large  injection  rates 

•  For  relatively  small  injection  rates  the  shaping  produces  a 
significant  stabilization  effect  in  the  mid-  and  far-field  relaxation 
regions 

•  The  suction-blowing  of  zero  mass  addition  destabilizes  the  flow 
in  the  whole  relaxation  region 
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Backup 
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Wall  heat  flux  (W/mA2) 


Stability  analysis  (cont’d) 


Run  2540 


Run  2540 


14 

12 

10 

8 

6 

4 

2 

0 


With  small  correction  of  y,  N-factors  predicted  by  the  perfect- 
gas  model  are  close  to  N-factors  predicted  by  STABL* 


*Wagnild,  R.M.  et  al.  AIAA-2010-1244 
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N  Factor 


j  Mean  flow  for  baseline  configuration 

(injection  rate=13.5  g/s) 


axial  velocity 


Amoth, 


Injection 

region 


1.00E+00 
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Mean  flow  near 
injection 

baseline  configuration 
injection  rate=13.5  g/s 


axial  velocity 


temperature 


Injection  region 


0  P2  0  13  Q14  D  ll 


OP?  a  IQ  0 19 
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It  is  not  easy  to  convert  temporal  growth  rates  to 

spatial  ones 


AMOTH, 


Gaster  transformation 


a  —  — 


V 

9 

V  —  du)  /da 

does  not  work  here 


jtOFSc^ 
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Spatial  stability  analysis  in  the  relaxation  region 

(iniection  rate^3.5  g/s,j^0.3) 


•Mode  0  (Mack  second  mode)  is  most  unstable 
•Its  instability  is  observed  at  low  phase  speeds 
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Injector  of  conical  shape 


/\MO>TH 


5-deg  half-angle  sharp  cone  with  the  injector  having 
the  slope  0=0C-A0 


A r  «  6 [tan  9  —  tan(0_  —  A0)]  «  bA9 
A  r  A  r 


A  b  = 


tan(0  A(9)  tail  9 
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Mean-flow  axial  velocity  for  conical  injectors 

(injection  rate  13.5  g/s) 


Baseline  (A0=O°) 


fli:  lOll  014  0t4  016 


Ott  OH  OX  Oil  022  OX  OX  OX  OX  027  OX  OX 


A0=1° 


13  014  016  ON 


OH  0t9  OX  Oil  OX  023  034  OX  OX  027  OX  OX 


A0=3° 


Oli  >0  13  014  OH  OH  >017  0 


019  OX  Oil  022  023  OX  OX  OX  027  OX  OX 


1  006+00 


-1  236-02 


1  ooe+oo 


0  006*00 


1.006*00 


Injection 

region 
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Shaping  of  A0=3°  with  injection  6.75  g/s 


Conical  injector  with  A6=3°  shaping 

•  Slightly  stabilizes  flow  in  the  near-field  relaxation  region 

•  Almost  zero  effect  in  the  far-field  relaxation  region 
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Shaping  of  A0=3°  with  injection  13.5  g/s 


Conical  injector  with  A6=3°  shaping 

•  Slightly  destabilizes  flow  in  the  near-field  relaxation  region 

•  Stabilizes  flow  in  the  far-field  relaxation  region 
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e  are  dealing  with  acoustic  instability 

(cont’d) 


•  Phase  speeds  of  unstable  modes  are  close  to  those  of  slow  acoustic  waves 

•  Resonant  interaction  can  enhance  receptivity  to  slow  free-stream  noise 

•  Instability  is  observed  at  low  frequencies  where  free-stream  noise  is  higher 
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i 


Maximal  growth  rates  in  the  relaxation  region 

with  13.5  g/s) 


0.24 


X 


•  The  most  unstable  mode  is 
Mack  second  mode  (mode  0) 

•  Unstable  x-region  decreases 
with  the  mode  number 


length  scale=8* 

o  —  —a  —  max  — or.fo;) 

max  fc.max  [  %  ^  ^  J 


Further  analysis  is  focused 
on  the  Mack  second  mode 
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max 


'gjjp , 

b>  Maximal  growth  rates  of  Mack  second 


0.0015 


0.0010- 


0.0005- 


0.0000 
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Run  2540 


length-scale  =  /  U 


x,  m 


In  the  relaxation  region  x>0.2  m 

•  High  maximal  growth  rates 

•  Low  frequencies 

•  Low  phase  speeds 
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Injection  affects  growth  rates  and 
frequencies  in  the  relaxation  region 


'£ 


Q) 

CD 


2 

O) 


T5:  Run  4520  Run  4520 


Maximal  growth  rates  are  increased 
Unstable  region  is 

•  narrowed  down  for  x<0.6 

•  widened  for  x>0.6 


Frequencies  are  decreased 
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Maximal  growth  rates  for  cylindrical  injector 

(various  injection  rates) 
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